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Abstract 



An experimental study of liquid drop impacts on a granular medium is proposed. 

Four fluids were used to vary physical properties: pure distilled water, water with 

glycerol at 2 concentrations 1 : 1 and 1 : 2 v/v and water with Tween 20 at the 

concentration of O.lg.l -1 . The drop free fall height was varied to obtain a Weber number 

(We) between 10 and 2000. Results showed that obtained crater morphologies highly 

depend on the impacting drop kinetic energy (Ek)- Different behaviours during the drop 

spreading, receding and absorption are highlighted as function of the fluids viscosity and 

surface tension. Experimental absorption times are also commented and compared with 

a simplified theoretical model. Drops maximal extensions and craters diameters were 
i I 

found to scale as We'-> and Et respectively. In both cases, found dependencies are 

smaller than those reported in literature: Wei for drop impacts on solid or granular 
i 

surfaces and E^ for spherical solid impacts on granular media. 
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1 Introduction 



The phenomenon of droplet impact on non-cohesive granular media is present in var ious in- 
dustr ial processes as food engineering, candy coated pills fabrication, starch wetting 



2008 



Vaclavik and Christi 



It is also of great importance in e nvironmental sciences: water erosion of s oil, pollution 

Ries et al.l . 12009] . . . However, 



of rivers and lakes, rainfall simulations jvan Diik et al.l. 12 002 



studies of this phenomenon are very scarce. [Holman et al. 



20021 studied the micro-metric 



drop impact on micro-metric grains of same sizes. Besides, authors used a very particular 
fluid, an aqueous polymer solution in this case. Others studied drops imbibition into powder 
beds with no impact . The drops were deposited on the granular surface with no kinetic en- 
ergy Hapgood et al.l . 12002] . Very recently, the problem has been investigated by Katsuragi . 



2010] where the author only focused on the granular medium deformation, in particular the 



resulting craters morphologies for different drops velocities and granular media (grains sized 
between 4 and 50/im). Moreover, he reported the existence of different craters types de- 
pending on the drop free fall height and that water drop maximal extensions D mr ,sr., that 



he as sumes equal to the crater diameter, scales as We 1 ^. However, Nefzaoui and Skurtys 



20091 



reported that, for water drops impact on coarser grains (average particle diameter 
81/im), D max ~ We 1 ^ and that the crater diameter D r i m is not equal to D max . 
On the other side, the drop impact phenome non onto thin liquid films or dry solid surfaces 
has been thoroughly studied for over a century Yarinl . 12006] . The drop behaviour after colli- 
sion with the surface is very complex since it depends on the physico-chemical characteristics 
of the drop and of the impact surface. The drop impact phenomenon can be divided in several 



sub-processes id e ntified as sp reading, receding, splashing and bouncing |Rioboo et al.l . 12002 



Wal et al.l . 120061 lYarinl . 12006] . In recent years more atte ntion has been shown to study the 



role o f the surface roughness on the drop impact process |Kannan and Sivakumarl . I2008L IXu 



2007] . In particular, it has been shown that the splashing of impacting drops is promoted by 
increasing the mean surface roughness. 

Spherical solid impact into granular media were investigated to some extent decades ago 



iNelson et alll2008[[Caballero et alll2007llWalsh et alll2003llZheng et aill2004ULohse et al. 



2004J. For impact energi es ranged between 10 2 and 0.6 J, the morphology of the craters 



was reported in detail by Walsh et all 12003] . In particular, the scaling of crater dimensions 
has been studied, and a power-law relationship between crater diameter D c and the energy 
of impact Ek can be derived in certain limits: D c ~ E*. The process of crater formation is 
complex since the resistance of the granular medium to the projectile penetration depends 
on the medium intrinsic properties, its packing mode and the applied force. Indeed, impact- 
ing sphere penetration dynamics and grain ele ction have b e en sh own to be very different 
whether the granular material is loose or dense |Lohse et al. . 2004 1. Anoth er type of exper- 



imen ts concerns a bead impact on a granular target made of similar beads |Beladjine et al. 



2007 



In this paper, an experimental study of drop impacts on a granular medium (glass beads) 
is proposed. Four different liquids were used in order to study the effect of viscosity and 
surface tension on the drop deformation for different impact velocities and try to determine 
the relevant parameters that governs the phenomenon. Different behaviours are highlighted 
depending on the receding magnitude and occurrence of splashing. The impact dynamics 
are also described through the spreading drop diameter temporal evolution obtained by high 
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speed video. Maximal spreading drop d i amete rs are compared to the case of drop impacts on 



solid surfaces |Clanet et al.l . 120041 iReinl . 11996] . Crater morphologies obtained by varying the 



drop kinetic energy at impact over more than three orders of magnitude are also presented 
and discussed. A power law for craters diameters dependence on the drop kinetic energy at 
impact is proposed and compared to solid projectiles impacts on granular media results. 



2 Material and methods 
2.1 Fluids and drop generation 

In our experiments, liquid drops were generated by a precision flat tipped syringe needle 
(Sigma-Aldrich, USA, St. Louis) connected to a digitally controlled syringe pump (Model 
1000, New Era Pump System Inc., Farmingdale, NY, USA). Liquid flow rate was sufficiently 
low to get a nil drop initial velocity (0.05 ml.min -1 ). Four fluids were used: pure distilled 
water, water with glycerol at 2 concentrations 1:1 and 1:2 v/v and water with Tween 20 at the 
concentration of O.lg.l" 1 . Their physical properties were measured and are presented in Table 
[TJ Three flat tipped needles (gauge 16, 18, 20) were used to obtain various drop diameters 
D , between 2.88 and 3.65 mm. Drops are assumed with some assurance to maintain a 
spherical shape through the free fall since their radius were lower than the capillary length k 
defined by . f^- |Landau and Lifshitz . 1959 1. This was verified on the captured images for the 



biggest drops. Different impact velocities V$ within the range of 0.4 to 5 m.s" 1 were obtained 
changing the drop release height h up to 1300mm. For each impact, Vq was measured from 
captured images. The kinetic energy of the drop at impact, given by Er = ^pD^gh was 
varied between 1 x 10 -6 and 3 x 10~ 4 J. Surface energy of the free drop surface, Eg = i\D\a 
was varied between 1.8 x 10~ 6 and 3 x 10~ 6 J. 



2.2 Granular medium 

Target surface was composed of a circular container which was 78 mm in diameter and 14 
mm deep. The container was filled with glass beads (G8893, Sigma-Aldrich, USA, St. Louis) 
which had an effective average particle diameter D g ~ 81/im and a density p g about 2300 
kg.m~ 3 . Granular medium compaction had a great effect on the crater formation. It was thus 
important to have the same initial beads bed state in general, and compaction in particular, 
for each drop impact. For this purpose, the same protocol was respected for each sample 
preparation. Beads were gently rolled into the container using a glass funnel. The container 
was overfilled and the surface was levelled using a straight-edge. There is no doubt that this 
operation introduced local compaction non-homogeneities. Finally, samples were weighed 
and the same mass of beads was used for each. Packing fractions of about 0.60 ± 0.02 were 
deduced from the granular material measurements of the mass and volume of the bulk. 



2.3 Fluid- granular contact angle 



The contact angle of water over the granular media was measured through a capillary suction 



experiment as detailed by |Xue et al. . 2006 1. The method consists in observing the velocity 
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Fluid 


p (kg.m 3 ) 


H (mPa.s) 


7 (mJ.m 2 ) 


Do (mm) 


Oh xlO 3 


Water 


996 


1 


72.2 2i 


58 - 3.25 - 3.65 


1.95 - 2.2 


Water-glycerol 1:1 


1121 


6 


66.9 


3.47 


11.8 


Water- glycerol 1:2 


1162 


19 


65.1 


3.39 


37.5 


Water-tween 20 (O.lg.l" 1 ) 


996 


1 


36 


3 


3 



Table 1: Measured physical properties of fluids, the drop diameter (D ) and the Ohnesorge 
number (Oh = ^J^ Do ) 



of a fluid rising in a granular column due to capillary forces. Taking into account hydro- 
static, viscous and capillary effects (only inertial effects are not considered here) leads to the 
following relation between the contact angle and the other physical parameters and variables 



XueetaL, 2006 



7 2 cos 2 (#) 



2pgp 



'3 

Z V 



3 z-2tv 

where 7, p, p, and 9 are the fluid surface tension, density, dynamic viscosity and the fluid- 
granular advancing contact angle respectively. Moreover, t and z are time and the rising fluid 

dz 

front position and v is the fluid front velocity given by v = — . The only unknown is 9 and 

at 

can easily be deduced from Eq. (CQ) since other parameters are measured experimentally. 



2.4 Granular medium pore diameter 



The granular medium mean po re diameter is needed to estimate the drop absorption times. 



According to |Xue et al. . 2006 1, Washburn equation with hydrostatic effects leads to the 



following expression of the pore average radius R p : 



1 



~fcos{6) 



t 

2rj 



pgz- 



(2) 



3 7 2 cos 2 (fl) 

This relation leads to consistent results when inertial effects are negligible, thus not at 
early times when t and z are too small. 



2.5 Images acquisition 

The impacting drop and the granular medium evolution, e.g. the drop spreading and receding, 
crater formation, were observed experimentally with a high-speed video Pulnix TM-6740GE 
(Pulnix, Inc., San Jose, CA , USA) with a pixel resolution of 640 x 480. It could capture 
200 frames per second (fps) in full frame and 1250fps with a 224 x 160 reduced matrix 
resolution. The camera was mounted on a boom stand which provided easy vertical or 
horizontal movement and stable support for the camera. A zoom video lens (18 — 108 mm 
f/2.5D, Edmund Optics, NJ, USA) was mounted on the video camera. Resulting frames were 
processed using ImageJ software (National Institutes of Health, USA). 
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3 Results and discussion 



A mathematical formulation of the impact leads to consider the equation of continuity and 
the momentum equations both in the radial and in the axial directions. Three dimensionless 
numbers can be considered after the adimensionalization procedure: Reynolds (Re = p ° D ° ), 

Weber (We = pD ^° ), Froude (Ft = In all tests, gravitational effect was negligible 

compared to inertia (Fr^ 1 <C 1). Thus, only the Weber number, which compares the inertial 
forces to surface tension and the Reynolds number (Re), ratio of inertial to viscous forces 
were retained as parameters of the problem. Thus, the dynamics at impact was driven 
by an interplay between the kinetic energy, viscosity and the droplet surface tension. The 
adimensional number -7= can also be used to take simultaneously into account the three 



parameters ( I ^ ertm - v ^ scoslt y \ Fedorchenko et al. 12005 1. The three typical phases observed 

1 v burjacel ension ' ' 1 J c r 



after the drop impact on the granular medium, i. e. the spreading, receding and absorption, 
are detailed in the next paragraphs. 



3.1 Drop spreading dynamics 

In Fig. [U the temporal evolution of the spreading factor (3 = D/D after impact is presented 
for high and low Weber numbers and different considered fluids. In a usual manner for 
inertia governed impacts, time t is made nondimensional using the impact velocity Vq and 
the initial spherical drop diameter D . In all cases, three phases can be observed: drop 
spreading > 0), receding < 0) and absorption. Spreading was the first stage of drop 
deformation after impact. The droplet initially flattens and spreads out horizontally into a 
pancake shape. Even if the description of the drop deformation before t = 0.8ms was beyond 
our experimental capacities, it is tantalizing to say that in the first instants, the drop was not 
deformed and just penetrates into the granular medium like a spherical ball for early times 
as shown in Fig. [3] (for t « 0). Indeed, in almost all first impact pictures, non-deformed 
spherical drops partially penetrating the granular medium were observed. This is not the 
case for a drop impacting on a solid surfac e where the d r op de formation starts at impact 
time and its diameter evolution scales as |Rioboo et al. . 2002|. 



Receding was the second stage of drop deformation after impact (see Fig. [I]). In the case 
of water and Tween drops, the receding behaviour was considerably different at intermediate 
(and at low We all the more) and high We. For We ~ 150, a "total receding" was observed, 
the drop final spreading factor j3j is relatively small: Df ~ IADq (water) and Df ~ 1.7Dq 
(tween). Besides, when it reaches its final diameter, the fluid recovers a drop-like shape (see 
Fig. [3]). Actually, at low and intermediate We, the drop spherical shape recovery let think 
that receding was governed by capillary forces. Whilst for higher kinetic energy, We ~ 750 
(water) and We ~ 590 (tween), final drop spreading factor was closer to (3 max : Df ~ 3.5D 
(water) and -D/ — 2.9D (tween). In contrast with the previous situation, we can talk about 
"partial receding" (see Fig. HI for t > 8 ms). In this case, the kinetic energy of the spreading 
fluid dominated the drop surface energy therefore the fluid groups into fingers along the 
outer rim of the droplet which had sufficient energy to overcome the surface tension and 
pinch off into smaller droplets. It is interesting to notice that a lower surface tension permits 
to increases surface area and hence D max . In fact, for an impact at We = 590, a Tween 
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drop reached the same maximal extension than a water drop impacting at We around 20% 
higher (We = 590) (see Fig. [1]). On the other hand, for highly viscous fluids (water-glycerol 
mixtures) with a surface tension close to water, the impacting drop receded to /3 < 1. Unlike 
water and Tween drops, no significative difference on the final drop diameter (after receding) 
was observed between high and low We impacts (see Fig. [1]). However, D max increases with 
We. Besides, oscillations of the spreading drop diameter were observed for low We. 

For the different considered fluids, splashing occurs at very high kinetic energy. To delin- 
eate the boundary between splashing and non-splashing regions a power-law relation based 
on Oh and Re numbers can be established: Oh x Re n = b where b is a real coefficient. Figure 
[2] shows Re and Oh values corresponding to the appearance of splashing. The boundary 
between splash-no splash then follows the relation: Oh x -Be 0,8 4 = 70 . This correlation can 
be compared with those reported by |Kannan and Sivakumar . 2008 | for a drop impinging 

0.85, OhxRe 117 = 63 



,0.61 



|Kannan and Sivakumar 
upon a dry solid surface or upon a thin fluid film where Oh x Re' 
respectively. As bfu m ~ b granu i ar , both exponents n can be compared. To splash, the same 
drop impacting on a granular medium needs higher kinetic energy than those impacting upon 
a thin fluid film since n f u m > n granu i ar . 









i i 

■ Water(We=130) 
-•- Water (We= 750) 

♦ Water-Tween20 (0.1 g.H) (We=165) 
— »— Water-Tween20 (0.1 g.H) (We=590) 
□ Water-Glyceroll:l (We=170) 

■ Water-Glyceral 1:2 (We=140) 
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Figure 1: Temporal evolution of the spreading factor (3 
and fluids. 



D j Dq for different Weber numbers 



3.2 Maximum drop deformation 



In Fig. [61 the maximal spreading factor (3 max = D max /DQ of drops impacting with different 
kinetic energies, is plotted versus We. All data were found to collapse in a 1/5 power law: 
Pmax We T o in spite of the surface tension and viscosity differences. A law of the type 

^ 2 1 — 2 

fimax oc We$ leads to (3 max cx Res since Wes oc D 5 Re$ . Our resu lts do not concord with 
the models describing the drop impact on a solid surface reported b y |Clanet et al. , 2004 1 and 
|Rein . 1996 1, (3 max oc We^ and (3 max oc Re^, respectively. Recently, Katsuragi . 2010 1 reported 
that after the impact of a free-falling water drop onto a granular layer (where the grain size 
is varied between 4 and 50/im), (3 max scales as We*. In the case of a unique drop and even 
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No splash 



• Water {D=2.88mm 
O Water [D=3.25mm 

• Water (D=3.65mm 
□ Water-Glycerol 1:1 
■ Water-Glycerol 1:2 

• Water-Tween 20 [0.1 g.l 
Oh=70xRe- 116 (R 2 =0.99) 



— 

le+03 
Re 



Splash 




Figure 2: Reynolds plotted versus Ohnesorge number showing the splash and no splash 
boundary. 

for drops of different diameters among values considered in this work ((D01/-D02) 5 G [1, 1-1]), 

1 2 . 

we can write We$ oc -Res. The insert of Fig. [6] confirms our relation between the maximum 
spreading factor /3 max and the Reynolds number Re for the four fluids with a high coefficient 
of determination (R 2 > 0.95). 



3.3 Drop absorption 

Figures [3] and H] show time sequences of a 2.88mm and 3.65mm diameter water drop impact 
for We ~ 100 and We ~ 650, respectively. Absorption is considered to be achieved when no 
visible difference between two and more consecutive frames can be observed. An absorption 
time of 128ms was obtained in the first case and 90ms in the second. It was surprising 
at first sight to find that the biggest drop took less time to be absorbed (an average bulk 
flow almost three times higher). After observing thoroughly the detailed frames sequences 
(which can not be shown here due to the high frames number), it has been found out that 

the mechanisms were not the same in the two cases. In the first case in particular, and for 

We 

all impacts at ^ 1 in general, the majority of drop absorption takes place after the 



Re 

drop spreading and receding and while the fluid is at rest. A fluid in a drop-like shape is 

then absorbed. Absorption and receding are quite independent and the latter is governed 

by capillary forces. In the second case, absorption is observed to start relatively early and 

We 

before the drop starts receding. This was observed for almost all impacts at , ^ 1 



Re 

where a flat film is absorbed instead of a spherical drop. In fact, the higher the impact 
kinetic energy, the bigger the drop maximal extension and the thinner the spreading film. 
Besides, the fluid film absorption does not start at the periphery, where the film thickness 
may be thought to be the lowest. Actu ally, the fluid film thickness decreas e s from the center 



to periphery at first spreading stages Clanet et al. . 2004 | |Roisman et al. . 2002 1. However 



when the spreading lamella reaches its maximal extension, a blob appears at the rim. While 
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Figure 3: Water drop impact images time sequence (Dq 
t m Oms, the top and front views are showed. 



2.88 mm; We ~ 150). Note: for 



there is no global fluid motion = in Fig. [T]), a local flow from the center to the pe r iphery 



still exists, feeding the blob and thinning the fluid film at its center |Roisman et al.l . 12002 



Fluid central sucking induces lamella receding. Therefore absorption is not the last stage 
of drop deformation. Absorption and receding are tightly coupled and we talk about an 
absorption induced receding. Table [2] reports measured absorption times for different fluids 
and for different release heights. A few interesting remarks can be made. 
First, these figures confirm that absorption times are systematically (except for very viscous 
fluids) higher for low release heights as explained above. Second, this difference of absorption 
times for different release heights decreases when the fluid viscosity increases. Finally, it 
vanishes for very viscous fluids (Fig. [5]). In fact, for glycerol-water 1:2 (19 times more 
viscous than water and water- tween mixture), no difference in absorption times is observed 
for different release heights while it can vary from one to 5 for water or water-tween mixture. 



It can be interesting to calculate theoretically absorption times and to compare them 
wit h the experimental res ults. The method presented and confirm ed in some particular cases 
by |Hapgood et al. . 2002 was adopted. In |Hapgood et al. . 2002 1 experiments, the drop was 
deposited with a syringe on the surface of the granular medium with no kinetic energy. With 
this method, the drop absorption was assumed to take place over a constant drop-granular 
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Figure 4: Water drop impact images time sequence (Dq = 3.65 mm; We ~ 650). 
medium contact surface (constant drawing area). Absorption time is then given by: 



Tcda 



VI 
1.35 d 



V 



e 2 R p >ycos(9 d ) 



(3) 



where Vd is the initial drop volume, e is the porous medium porosity, R p is the pore radius, 7, 
77 and 9d are the fluid surface tension, viscosity and the fluid-granular medium contact angle 
respectively. 

In our case, calculated times corresponded to experimental times in very few cases (See 
Table |2]), for water drops falling from very low heights for example (experimental conditions 
very close to those of Hapgood) and no concluding results were obtained for water-tween and 
water-glycerol mixtures. The calculated times depended very highly on contact angles which 
were quite difficult to determine with a good precision. Besides, using this model implies 
the implicit assumption that no absorption occurs during spreading and receding which is 
obviously not true. 
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— &8 tm 




f > 1248 ms 



Figure 5: Water-glycerol mixture (bulk ratio 1:1) drop impact images time sequence (Do 
3.47 mm; We ~ 570). 



3.4 Crater morphologies 

In Figure [71 for different Weber numbers, examples of crater morphologies observed by im- 
pacting into a container of glass beads a 3.65 mm diameter water drop are shown. The 
granular medium keeps a memory of the fluid shape at its maximal extension, thanks to its 
deformability and wettability. Similar craters were also observed with the different fluids. 
At low We, e.g. We ~ 30 (Fig. 0-A), a crater with in its centre a circular agglomeration of 
glass beads was formed. A rotationally symmetric rim was raised above the original target 
surface level. Glass beads driven by the receding drop formed a central agglomeration where 
the surface energy is sufficiently strong to join glass beads by capillary cohesion forces. This 
crater morphology was observed for We < 100. As We was increased the first irregularities 
appeared inside the crater and the agglomeration in the centre was less circular (Fig. [7]-B). 
These irregularities became more marked at higher We (Fig. [7J-C-D). They can be associ- 
ated with the fingering instability at the rim of the ex panding water drop like reported in th e 



ated witn tne nngermg instability at tne rim ol tne ex panding water drop like reported m tn 
previous studies of the water drop on a solid surface |AllenlTl975 . Mehdizadeh et al. . 2004 



This crater morphology occur red in the range 100 < We < 300. This kind of crater was not 
reported by Katsuragil . |2010| . In its work, finer grains were used and lower We numbers 
were considered which may explain these differences. At higher kinetic energy, a layer of wet 
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Figure 6: The maximal spreading factor /3 max = D max /D as function of the Weber number 
(We) for different drop diameters and fluids. Insert: The maximal spreading factor (3 max = 
Dmax/Do as function of the Reynolds number, Re. 



Fluid 


D(mm) 


We 


t(s) 


r(s) 


Water 


2.88 


150 


0.12 


0.128 


Water 


3.65 


750 


0.09 


0.09 


Wat er-t ween 


3 


164 


0.27 


1.25 


Wat er-t ween 


3 


549 


0.05 


1.25 


Water-Glycerol 1:1 


3.47 


114 


1.68 


6.3 


Water-Glycerol 1:1 


3.47 


570 


1.06 


6.3 


Water-Glycerol 1:2 


3.39 


119 


2.29 


14.1 


Water-Glycerol 1:2 


3.39 


546 


2.43 


14.1 



Table 2: Measured (t) and calculated (r) absorption times for different drops of different 
sizes and fluids and falling from different heights. 



glass beads appears instead of the central agglomeration previously described (Fig. 0-E). It 
is due to the water drops splashing. 

In Figure [HI the crater rim diameter D rim , measured at the top of the crater rim, is 
plotted against kinetic energy of drops at impact. As expected, for different drop diameters 
and fluids, the measured diameter D r i m increased with increasing Ek- Moreover, the data are 
well described by a | power law over more than three orders of magnitude in Ex- However, 
this law is slightly different of those reported for the impact of a solid projectile where D rim 

scales as Et IWalsh et al.. 



(Katsuragil . 



20031 lUehara et al.l . 12003] and the water drop impact reported by 



2010 
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(E) We a 380 

Figure 7: Crater morphologies observed by dropping a 3.65 mm diameter water drop into a 
container of glass beads for different Weber numbers. 




Figure 8: Crater rim diameter (D r i m ) as function of the kinetic energy (Ek) for different 
drop diameters and fluids. 
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4 Conclusion 



An experimental study of liquid drop impacts on a granular medium has been presented 
for Weber numbers (We) in the range 10 — 2000. Four liquids were used to vary either 
the viscosity (between 1 and 19mPa.s) or the surface tension (between 36 and 72mJ.m~ 2 ). 
Different crater morphologies have been observed depending on the Weber number. The 
crater rim diameter was found to scale as a 1/5 power-law of the kinetic energy of drops 

at impact Ek- This law was slightly different of those reported for the impact of a solid 

i 

projectile where D r i m scales as E^. 

The temporal evolution of the spreading factor (/3 = D/Dq) was presented. The effects 
of viscosity, surface tension and inertia on the spreading, receding and absorption were com- 
mented. A pow er-law (3 max = D max /D ~ We^ smaller than those recently reported by 



Katsuragil |2010| has been deduced for the different fluids. 

The boundary between splash-no splash was found to follow the relation: OhxRe 84 " = 70. 
In particular, to splash, the same drop impacting on a granular medium needs higher kinetic 
energy than those impacting upon a thin fluid film. 

Measured absorption times depended highly on We and the fluid vis cosity. Calculated 



times according to the simplified model reported by lHapgood et al.l |2002| were not in agree- 



ment with experiments. A better prediction of absorption times must be achieved taking 
into account absorption during spreading and receding. 

Finally, more imagery acquisition is necessary to understand the fluid granular interactions 
(crater formation mechanisms) and to develop an energetic model that details the initial drop 
kinetic energy dissipation during the drop deformation (viscous dissipation, . . . ) and crater 
formation (granular dissipation). 



Acknowledgements 

The authors thank F. Melo for discussions and help, the CONICYT sponsor of the Bicente- 
nario Grant (PBCT PSD-62) and Dicyt (project 08107150). 



14 



References 



R.F. Allen. The role of surface tension in splashing. J. Colloid Interface Sci., 51:350, 1975. 

D. Beladjine, M. Amni, A. Valance, and L. Oger. Collision process between an incident bead 
and a three-dimensional granular packing. Phys. Rev. E, 75:061305, 2007. 

G. Caballero, R. Bergmann, D. Van der Meer, A. Prosperetti, and D. Lohse. Role of air in 
granular jet formation. Phys. Rev. Lett., 99:018001, 2007. 

C. Clanet, C. Beguin, D. Richard, and D. Quere. Maximal deformation of an impacting drop. 
J. Fluid Mech., 517:199, 2004. 

A. I. Fedorchenko, A.B. Wang, and Y.H. Wang. Effect of capillary and viscous forces on 
spreading of a liquid drop impinging on a solid surface. Phys. Fluids, 17(9):093104, 2005. 

K.P. Hapgood, J.D. Litster, S.R. Biggs, and T. Howes. Drop penetration into porous powder 
beds. J. Colloid Interface Sci., 253(2) :353, 2002. 

K.R. Holman, M.J. Cima, S.A. Uhland, and E. Sachs. Spreading and infiltration of inkjet- 
printed polymer solution droplets on a porous substrate. J. Colloid Interface Sci., 249(2): 
432, 2002. 

R. Kannan and D. Sivakumar. Impact of liquid drops on a rough surface comprising mi- 
crogrooves. Exp. Fluids, 44:927, 2008. 

H. Katsuragi. Morphology scaling of drop impact onto a granular layer. Phys. Rev. Lett., 
104:218001, 2010. 

L.D. Landau and E.M. Lifshitz. Fluid Mechanics, Course of Theoretical Physics Vol. 6. 
Addison- Wesley, Reading, 1959. 

D. Lohse, R. Bergmann, R. Mikkelsen, C. Zeilstra, D. Van der Meer, M. Versluis, K. Van der 
Weele, M. Van der Hoef, and H. kuipers kuipers. Impact on soft sand: Void collapse and 
jet formation. Phys. Rev. Lett., 93:198003, 2004. 

N.Z. Mehdizadeh, S. Chandra, and J. Mostaghimi. Formation of fingers around the edges of 
a drop hitting a metal plate with high velocity. J. Fluid Mech., 510:353, 2004. 

E. Nefzaoui and O. Skurtys. Water drop impact on granular medium. In Book of Abstracts. 
Southern Workshop on Granular Materials - SWGM09, Nov. 30 to Decs 4, 2009 - Vina 
del Mar, Chile, 2009. 

E.L. Nelson, H. Katsuragi, P. Mayor, and D.J. Durian. Projectile interactions in granular 
impact cratering. Phys. Rev. Lett., 101:068001, 2008. 

M. Rein. The transitional regime between coalescing and splashing drops. J. Fluid Mech., 
306:145, 1996. 



15 



J.B. Ries, M. Seeger, T. Iserloh, S. Wistorf, and W. Fister. Calibration of simulated rainfall 
characteristics for the study of soil erosion on agricultural land. Soil and Tillage Research, 
106(1):109, 2009. 

R. Rioboo, M. Marengo, and C. Tropea. Time evolution of liquid drop impact onto solid, 
dry surfaces. Exp. Fluids, 33:112, 2002. 

Ilia V. Roisman, Romain Rioboo, and Cameron Tropea. Normal impact of a liquid drop on 
a dry surface: model for spreading and receding. Proc. R. Soc. London, Ser. A, 458(2022): 
1411, 2002. 

J.S. Uehara, M.A. Ambroso, R.P. Ojha, and D.R. Durian. Low-speed impact craters in loose 
granular media. Phys. Rev. Lett, 90:194301, 2003. 

V.A. Vaclavik and E.W. Christian. Essentials of Food Science. Springer, third edition, 2008. 

A. I. J. M. van Dijk, L. A. Bruijnzeel, and C. J. Rosewell. Rainfall intensity-kinetic energy 
relationships: a critical literature appraisal. J. Hydrol. (Amsterdam, Neth.), 261(1-4):1, 
2002. 

R.L. Vander Wal, CM. Berger, and S.D. Mozes. The splash / non-splash boundary upon a 
dry surface and thin fluid film. Exp. Fluids, 40:53, 2006. 

A.M. Walsh, K.E. Holloway, P. Habdas, and J.R. de Bruyn. Morphology and scaling of 
impact craters in granular media. Phys. Rev. Lett., 91:104301, 2003. 

L. Xu. Liquid drop splashing on smooth, rough, and textured surfaces. Phys. Rev. E, 75: 
056316, 2007. 

H.T. Xue, Z.N. Fang, Y. Yang, J. P. Huang, and L.W. Zhou. Contact angle determined 
by spontaneous dynamic capillary rises with hydrostatic effects: Experiment and theory. 
Chem. Phys. Lett, 432:326, 2006. 

A.L. Yarin. Drop impact dynamics: Splashing, spreading, receding, bouncing .... Annu. 
Rev. Fluid Mech., 38:159, 2006. 

X.J. Zheng, Z.T. Wang, and Z.G. Qiu. Impact craters in loose granular media. Eur. Phys. 
J. E: Soft Matter Biol. Phys., 13:321, 2004. 



16 



